Draft version April 24, 2013 

Preprint typeset using L^T^X style cmulatcapj v. 12/01/06 



REDSHIFT 6.4 HOST GALAXIES OF 10 s SOLAR MASS BLACK HOLES: 
LOW STAR FORMATION RATE AND DYNAMICAL MASS 

Chris J. Willott 

Herzberg Institute of Astrophysics, National Research Council, 5071 West Saanich Rd, Victoria, BC V9E 2E7, Canada 



o 

< 

(N 
(N 

o 

u 

6 

CO 

> 
00 

in 

T—i 

O 

m 

> 

X 



Alain Omont and Jacqueline Bergeron 

UPMC Univ Paris 06 and CNRS, UMR7095, Institut d'Astrophysique de Paris, F-75014, Paris, France 

Draft version April 24, 2013 

ABSTRACT 

We present ALMA observations of rest-frame far-infrared continuum and [C n] line emission in two 
z = 6.4 quasars with black hole masses of « 10 8 M Q . CFHQS J0210-0456 is detected in the continuum 
with a 1.2 mm flux of 120 ± 35 /zJy, whereas CFHQS J2329-0301 is undetected at a similar noise level. 
J2329-0301 has a star formation rate limit of < 40M Q yr _1 , considerably below the typical value at 
all redshifts for this bolometric luminosity. By comparison with hydro simulations, we speculate that 
this quasar is observed at a relatively rare phase where quasar feedback has effectively shut down star 
formation in the host galaxy. [Cn] emission is also detected only in J0210-0456. The ratio of [Cn] to 
far-infrared luminosity is similar to that of low redshift galaxies of comparable luminosity, suggesting 
the previous finding of an offset in the relationships between this ratio and far-infrared luminosity 
at low- and high-redshift may be partially due to a selection effect due to the limited sensitivity of 
previous continuum data. The [C II] line of J0210-0456 is relatively narrow (FWHM = 189±18kms _1 ), 
indicating a dynamical mass substantially lower than expected from the local black hole - velocity 
dispersion correlation. The [C n] line is marginally resolved at Of.' 7 resolution with the blue and red 
wings spatially offset by CK'5 (3 kpc) and a smooth velocity gradient of 100 km s -1 across a scale of 
6 kpc, possibly due to rotation of a galaxy- wide disk. These observations are consistent with the idea 
that stellar mass growth lags black hole accretion for quasars at this epoch with respect to more recent 
times. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift — quasars: 
general 



1. INTRODUCTION 

The peak of global star formation occurred about 10 
billion years ago at a redshift of z ~ 2 (Reddy & Stei- 
del 2009). The rise in star formation at earlier times is 
studied by tracing the space density and properties of 
higher redshift galaxies. Such galaxies can be selected 
in the rest-frame ultraviolet as Lyman break dropouts 
(Bouwens et al. 2006), in the rest-frame infrared as line 
or continuum sources (Carilli & Walter 2013), via black 
hole accretion activity as active galactic nuclei (AGN) 
or quasars (Fan et al. 2006), or via stellar explosions 
such as gamma ray bursts (Tanvir et al. 2012). These 
methods are complementary in that they are sensitive to 
galaxies with a range of mass, star formation rate, dust 
formation rate and black hole accretion rate, to allow a 
broad view of early galaxy evolution. 

Thanks to the high fraction of stellar radiation re- 
radiated in the infrared by interstellar dust and gas and 
the negative k-correction (Blain & Longair 1993), high 
redshift galaxies can be well studied by millimeter obser- 
vations. Continuum observations are sensitive to the star 
formation rate as radiation from young, hot stars is re- 
radiated by dust. Molecular lines such as CO probe the 
molecular gas in dense star-forming regions. Atomic lines 
such as the fine-structure line of singly-ionized carbon, 
[C n] , probe the interstellar medium and the outer parts 
of star-forming regions. It has been recognized that the 
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[C Ii] line will likely become the most useful line for study- 
ing very high redshift galaxies with the Atacama Large 
Millimeter Array (ALMA; Walter & Carilli 2008). In- 
deed early ALMA observations already show detections 
of [C n] in normal star-forming galaxies at z > 4 (Wagg 
et al. 2012; Carilli et al. 2013). [Cn] has also been 
detected in a few z > 6 quasar host galaxies (Maiolino 
et al. 2005; Walter et al. 2009; Venemans et al. 2012; 
Wang et al. 2013). 

One of most puzzling aspects of galaxy evolution is 
the tight correlation between black hole mass and galaxy 
properties such as bulge stellar mass and velocity disper- 
sion (Magorrian et al. 1998; Ferrarese & Merritt 2000). 
This correlation suggests a physical connection between 
black hole accretion on sub-parsec (pc) scales and galaxy- 
wide star formation and gas accretion on kiloparsec (kpc) 
scales. The most likely explanation is quasar feedback, 
but the details of how this operates as a function of 
cosmic time are still to be determined (Cattaneo et al. 
2009). Observationally, this topic can be studied by mea- 
suring the global growth rates of black holes and galaxies 
and the ratio of black hole to stellar mass in high redshift 
galaxies. The black hole masses of quasars can be mea- 
sured from the dynamics of gas in the broad line region 
which is only ~ 1 pc from the black hole (Wandel 1999). 

As well as studying the physical properties of star- 
forming gas in high-redshift galaxies, millimeter intcr- 
ferometry can be used to probe the gas dynamics. A 
particularly useful application is that of quasar host 
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galaxy dynamical masses to measure the ratio of black 
hole to galaxy mass in the early universe (Walter et al. 
2004). Observations in the rest-frame UV or optical are 
hampered by the overwhelming brightness of the quasar 
point-source (e.g. Mechtley et al. 2012). Wang et al. 
(2010) showed that CO line widths of the most optically 
luminous z 6 quasars indicate ratios of black hole to 
galaxy mass a factor of 10 on average greater than found 
in the local universe. This could signal that black holes 
grow much more rapidly than their host galaxies within 
the first billion years. However, there is a selection bias 
detailed in depth by Lauer et al. (2007) which suggests 
that the most luminous quasars will be biased to high 
black hole mass due to scatter in the correlation. In 
order to check whether this bias affects the conclusions 
of Wang et al. (2010) it is important to determine the 
black hole to galaxy mass relationship for the more com- 
mon, lower luminosity quasars at z sa 6 such as those 
identified in the Canada-France High-z Quasar Survey 
(CFHQS; Willott et al. 2010a). 

An alternative approach to studying the co-evolution 
of galaxies and black holes is to determine the star for- 
mation rate of active galaxies over cosmic time (Carilli et 
al. 2001; Omont et al. 2003; Priddey et al. 2003; Wang 
et al. 2008, 2011a; Lutz et al. 2010; Serjeant et al. 2010; 
Bonfield et al. 2011; Omont et al. 2013). These stud- 
ies suggest positive evolution in the star formation rate 
(SFR) at a fixed quasar luminosity from z = to z rj 2 
and approximately constant or a mild decline at higher 
redshift. However, in these studies, most high-redshift 
quasars are undetected and conclusions have to be based 
on statistical detections of stacked sub-samples. ALMA 
now provides an opportunity to measure the SFR at least 
an order of magnitude fainter than previous observations 
and revolutionize our understanding of the co-evolution 
of black holes and their host galaxies. 

In this paper we present an ALMA study of the in- 
terstellar medium of the host galaxies of two z = 6.4 
CFHQS quasars. These quasars were selected for study 
based on their high redshift (they are two of the four 
highest redshift published quasars) , faint absolute magni- 
tude (M1450 > —25) and black hole masses of ~ 10 s M Q . 
Section 2 details the new observations. The results are 
presented in Section 3. Section 4 contains a discus- 
sion of the results. Cosmological parameters of Hq = 
70 km s _1 Mpc -1 , fl M = 0.27 and Q A = 0.73 (Komatsu 
et al. 2011) are assumed throughout. 

2. OBSERVATIONS 

CFHQS J021013-045620 (hereafter J0210-0456) and 
CFHQS J232908-030158 (hereafter J2329-0301) were ob- 
served with ALMA between June and August 2012 in 
Early Science project 2011.0.00243.S. The number of 
12 m diameter antennae in use ranged from 17 to 24 with 
a typical longest baseline of 400 m. Observations of the 
science targets were interleaved with nearby phase cal- 
ibrators, J0217+017 and J2323-032. Uranus was used 
as the amplitude calibrator and 3C446 as the bandpass 
calibrator. Total on-source integration times were 8000 s 
for J0210-0456 and 8500 s for J2329-0301. 

The band 6 (1.3 mm) receivers were set up so that 
one of the four basebands (each of width 1.875 GHz) was 
centred on the expected location of the rcdshifted [C n] 
transition (^ rest =1900.5369 GHz). The redshifts adopted 
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Fig. 1. — ALMA 1.2 mm continuum images generated from the 
three line-free basebands for each of the two quasar fields. The 
greyscale ranges from — 3cr (black) to +3cr (white) where a = 
35, 30 fiJy beam" 1 for J0210-0456 (left) and J2329-0301 (right), re- 
spectively. Red circles show the expected source positions (circles 
do not indicate the positional uncertainty). There is a 3.4 a detec- 
tion for J0210-0456 and no detection for J2329-0301. The strongest 
source in the field of J2329-0301 is marked with a blue circle and 
labelled 'b'. It is co-incident with a blue galaxy in the optical 
images of Willott et al. (2007). 

were those of the low-ionization broad Mgn lines of the 
quasars measured by Willott et al. (2010b). Previous 
studies of high-redshift quasars have shown relatively 
small offsets (la dispersion 270 km s -1 ) between Mgn 
and the systemic redshift (Richards et al. 2002). The 
remaining three spectral windows were placed nearby to 
sample the 1.2 mm continuum. Each baseband is sam- 
pled by 120 channels of width 15.625 MHz (equivalent to 
wlSkms- 1 ). 

Data processing was performed by staff at the North 
American ALMA Regional Center using the CASA soft- 
ware package 1 . The three line-free spectral windows 
were combined to generate 1.2 mm continuum images. 
Both the continuum maps and spectral line datacubes 
were spatially sampled with O'.'l pixels. The synthesized 
beams are 0777 by 0752 for J0210-0456 and 0'.'73 by 0763 
for J2329-0301. The noise level reached in a 2 channel 
bin (31.25 MHz) is 0.22 mJy beam" 1 for J0210-0456 and 
0.23 mJy beam" 1 for J2329-0301. 

3. RESULTS 
3.1. Far-infrared luminosity 

The 1.2 mm continuum luminosity of z = 6.4 sources 
probes rest-frame 160 /im radiation, on the Rayleigh- 
Jeans tail side of the typical star-forming galaxy dust 
spectral energy distribution (SED; Lagache et al. 2005). 
This makes it an excellent proxy for the total far- 
infrared luminosity (Lfir; integrated luminosity over 
42.5 — 122.5 /im) which is a reliable tracer of the star 
formation rate due to dust heated by young stars. In 
the most ultraviolet-luminous quasars (such as those at 
z ~ 6 in the SDSS), there is often a substantial contri- 
bution to Lfir from dust heated by the AGN (Wang et 
al. 2008). The CFHQS quasars are an order of magni- 
tude less UV-luminous than SDSS quasars and therefore 
should have a correspondingly lower contribution from 
AGN- heated dust, allowing continuum observations to 
probe lower star-formation rates. 

The ALMA 1.2 mm continuum images generated from 
the three spectral windows that did not include the 
[C n] line were analyzed to determine their flux-densities. 
These images are shown in Figure [l] where the expected 

1 http://casa.nrao.edu 
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Fig. 2. — Far-infrared luminosity versus AGN bolometric luminosity for z ~ 6 quasars. The two CFHQS quasars observed with ALMA 
in this paper are shown with triangles. Previous detections from Wang et al. (2011a), mostly of SDSS quasars with a few quasars from 
other surveys, are shown as open circles. Previous detections of three CFHQS quasars are shown as filled circles (Omont et al. 2013). The 
squares show stacked averages from Wang et al. (2011a) and Omont et al. (2013). A complete sample of local PG quasars are shown with 
gray crosses (Hao et al. 2005). Colored lines show stacked averages in bins in redshift and LboI for quasars from the H- ATLAS survey and 
other far-IR/mm data (Serjeant et al. 2010). The dotted line is a fit to high-redshift (2 < z < 7) stack averages (Wang et al. 2011a). The 
dashed line converts LboI to black hole accretion rate and Lpi^ to star formation rate such that these grow in parallel to match the local 
A^BH/A^stcllar relationship (Tundo et al. 2007). 



locations of the quasars are identified by red circles. 
J0210-0456 is detected at 3.4 a with /i. 2mm = 120 ± 
35 /iJy. At this significance level it is not possible to de- 
termine whether the source is spatially resolved. J2329- 
0301 is undetected with no hint of positive flux at the 
quasar location. The most significant continuum source 
in the field (7" north of the quasar, labeled 'b') is identi- 
fied as a blue galaxy at much lower redshift in the optical 
imaging of Willott et al. (2007). 

This continuum flux-density was converted to a far- 
infrared luminosity assuming a typical SED for high- 
redshift star-forming galaxies. To make meaningful com- 
parison with previous results (in particular Wang et al. 
2008; 2011a, Omont et al. 2013) we adopt a greybody 
spectrum with dust temperature, = 47 K and emissiv- 
ity index, /? = 1.6. We note that our faint sources have 
much lower mm fluxes than the typical sources used to 
determine these parameters. If our sources instead have 
dust temperature closer to that of nearby luminous in- 
frared galaxies (LIRGs, 10 11 - 10 12 L ,T d w 33 K, U et 



al. 2012) then the values of Lfir would be 3x lower. 
For the remainder of this paper, uncertainties on Lfir 
(and inferred SFR) only include the flux measurement 
uncertainties, not that of the dust temperature. 

The far-IR luminosity of J0210-0456 is (2.60 ± 0.76) x 
lO n L . J2329-0301 is undetected with L Fm < 1.9 x 
10 11 L (3cr limit). We note the incredible sensitivity 
of these early ALMA observations that reach the lower 
end of the LIRG classification in the early universe at 
z = 6.4. 

We now consider the relationship between Lfir and 
the quasar bolometric luminosity LboI f° r high-redshift 
quasars. L^oi in this case is for the quasar component of 
the galaxy and assumes a typical bolometric correction 
from the rest-frame UV luminosity at 1450 A of a factor 
of 4.4 (Richards et al. 2006). Lb i does not include any 
excess FIR luminosity above that of the typical quasar. 
It is still a matter of debate as to how much of the typical 
quasar far-IR emission is due to dust heated by the AGN, 
compared to dust heated by a starburst (Haas et al. 2003; 
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Hao et al. 2005; Netzer et al. 2007; Lutz et al. 2010). 

Because most high-redshift quasars have not been de- 
tected in mm continuum with the sensitivity level of pre- 
vious studies, the relationship between Lfir and LboI 
has been based on stacking of sub-samples with differ- 
ent bolometric luminosity ranges. Wang et al. (2011a) 
showed that the stacks based on several samples at 
2 < z < 7 could be fit by the relationship Lfir °^ ^BoV 
Omont et al. (2013) found that the stacked average 
from 1.2 mm MAMBO observations of CFHQS z « 6 
quasars also lie on this relationship. The non-linear na- 
ture and significant scatter (for those detected so far) is 
interpreted in an evolutionary scenario where both star 
formation rate and black hole accretion are dependent 
upon dark matter halo mass, but with a lack of synchro- 
nization in the rates of these processes. The galaxies with 
the lowest Lfir at a given Lb i are expected to have a 
significant fraction of their Lfir due to quasar-heated 
dust (Netzer et al. 2007). 

Figure [2] shows previously published data for 
individually-detected z w 6 quasars and stacked aver- 
ages from Wang et al. (2011a) and Omont et al (2013). 
The dotted line is the relationship between Lfir and 
LboI found by Wang et al. (2011a) for stack averages 
of quasars at 2 < z < 7. Note that both the stacked 
averages and relationship from Wang et al. 2011a have 
been renormalized according to the bolometric correc- 
tion adopted here (see Omont et al. 2013 for more de- 
tails). Gray crosses are a complete sample of low-redshift 
(z < 0.5) optically-selected Palomar-Green (PG) quasars 
(Hao et al. 2005). Lfir for PG quasars has been es- 
timated as 2 x the luminosity at 60 jxm (Lawrence et 
al. 1989). Note that many of the highest luminosity 
(most distant) quasars in the PG sample are undetected 
at 60 /xm and only have upper limits on Lfir. As noted 
by Wang et al. (2011a), the z ~ 6 stacked averages lie 
close to the correlation exhibited by low-redshift quasars, 
indicating no enhancement in SFR at high redshift for a 
given quasar luminosity. 

Serjeant et al. (2010) used Herschel imaging of quasars 
in the H-ATLAS survey plus supplementary published 
IR and mm data to determine the average quasar far- 
infrared luminosity as a function of both redshift and 
quasar luminosity. Their data (for all bins containing 
10 or more quasars) is shown in Figure using an / 
band bolometric correction of 12.0 (Richards et al. 2006) 
and Lfir = 1-75 x the luminosity at 100 /im. These 
data show a similar correlation of the two luminosities as 
found for the PG and z ~ 6 samples. However, Serjeant 
et al. do find a positive correlation between Lfir and 
redshift up to z s» 3 that does not continue up to the 
z — 6 data of Wang et al. (2011a) and Omont et al. 
(2013). A positive correlation between Lfir and redshift 
up to z w 2 was also observed by Bonfield et al. (2011). 

The values for the two CFHQS quasars with ALMA 
data are also plotted on Figure[2j The detection of J0210- 
0456 shows its Lfir to be a factor of w 3 lower than the 
stacked average relationship and below the stacked aver- 
ages from H-ATLAS at all redshifts. The non-detection 
of J2329-0301 corresponds to Lfir at least a factor of 
10 lower than the stacked average from the full CFHQS 
sample (Omont et al. 2013) and substantially below the 
H-ATLAS averages. 

The far-infrared luminosity can be used to derive the 
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Fig. 3. — Continuum-subtracted [Cn] line emission for 
CFHQS J0210-0456 integrated over 15 channels containing the 
[Cn] line. Contours are at [—3,3,5,7] X a where <r = 
0.03 Jy kms - 1 beam -1 . 

SFR assuming the relation in Kennicutt (1998) with a 
Salpeter (1955) initial mass function (IMF). For J0210- 
0456, SFR= 48M (? yr- 1 and for J2329-0301, SFR< 
40M Q yr _1 (3 a limit). For both these objects, the as- 
sumption in deriving SFR is that there is no contribu- 
tion to Lfir from quasar-heated dust. These quasars lie 
close to the lower range of L F ir where it has been sug- 
gested that the majority of the cool dust is heated by the 
quasar (Netzer et al. 2007) . If this is the case then the 
SFR would be even lower. The very low SFR implied 
for the host galaxy of J2329-0301 is surprising given that 
it has a 2.5 x 10 s M Q black hole accreting at the Ed- 
dington rate (Willott et al. 2010b) and a very luminous, 
spatially-extended Lya halo (Goto et al. 2009; Willott 
et al. 2011). 

3.2. [Cll] luminosity 

The datacubes of the spectral windows containing the 
expected [C ii] emission lines for the two quasars were 
inspected for line emission. A line was easily detected for 
J0210-0456, but not for J2329-0301. Figure [3] shows the 
image of [Cn] emission for J0210-0456 integrated over 
the 15 channels (each of width 15.625 MHz) that show 
line emission. Continuum emission has been subtracted 
off using the continuum image of Figure [T] The source 
is elongated east- west, although note this is close to the 
major axis of the elongated beam. The spatial structure 
will be discussed further in the following section. No 
other [C Ii] emitters at the same redshift are seen in the 
field. 

The spectrum of J0210-0456 is plotted in Figure[I] The 
[C n] line is offset from the broad ultraviolet Mg n emis- 
sion line by 230 kms -1 . Note that the rms observational 
uncertainty on the Mgll redshift is 160 km s -1 , so the 
redshifts of [Cn] and Mgll are consistent. We take the 
redshift of Zfciii = 6.4323 ± 0.0005 to be the systemic 
redshift because it is measured to much higher accuracy 
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Fig. 4.— [Co] spectrum for CFHQS J0210-0456 overlaid with 
best fit Gaussian plus continuum model (blue). The red square 
with error bar is the continuum flux measured from the three line- 
free basebands. The arrow marked Mgll shows the redshift mea- 
sured from the quasar broad line region. 

than Mgll and is associated with star formation in the 
host galaxy rather than gas in the circum-quasar envi- 
ronment. A simple Gaussian plus constant fit was made 
to the observed spectrum. A Gaussian with FWHM 
= 189 ± 18 kins -1 provides a good fit to the line pro- 
file. The best-fit constant is positive showing a non-zero 
continuum level that is consistent with the continuum 
flux-density of 120 /iJy measured from the three line- free 
basebands in Section [3~T1 

The improved systemic redshift for J0210-0456 allows 
an improvement in the determination of the size of the 
highly-ionized near-zone. The size of the region ionized 
by the quasar depends upon several factors including the 
neutral hydrogen fraction of the intergalactic medium 
when the quasar first turned on (Madau & Rees 2000; 
Cen & Haiman 2000) and therefore can be used to probe 
cosmic reionization. Willott et al. (2010b) used the Mgn 
redshift of this quasar to determine a near-zone size of 1.7 
proper Mpc, which is lower than any other z ss 6 quasar 
except for linclcss or broad absorption line quasars (Fan 
et al. 2006; Carilli et al. 2010). Using the new redshift 
°f z [cii] — 6.4323 gives a near-zone size of 1.4 proper 
Mpc, even lower than previously calculated. Correcting 
for the known luminosity-dependence of R oc L 1 / 3 makes 
the size only slightly lower than the typical size for more 
luminous z > 6.1 quasars of 5 Mpc (Carilli et al. 2010). 

The spectrum of J2329-0301 is plotted in Figure [5) 
It can be seen that there is not strong evidence for a 
measurable [C n] emission line. There is weak positive 
flux at the Mgn redshift that may correspond to real 
emission, but it is very uncertain so we assume here a 
non-detection. 

The [C n] line flux of J0210-0456 was determined by 
integrating over the channels containing the line and sub- 
tracting off the continuum component. This line flux was 
then converted to a line luminosity at the measured red- 
shift. For J2329-0301 a 3(7 upper limit for the [Cn] flux 
and luminosity was determined by assuming a spatially 
unresolved Gaussian with FWHM = 300 km s -1 . This 
line width is somewhat broader than that observed for 
J0210-0456 but narrower than CO line widths for SDSS 
z s=! 6 quasars (Wang et al. 2010). Measurements of emis- 
sion line and continuum parameters derived from these 
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Fig. 5.— [Cn] spectrum for CFHQS J2329-0301. There is no 
convincing detection of the line for this quasar. The red square 
with downward arrow is the limit on the continuum flux from the 
three line-free basebands. 



TABLE 1 
Millimeter data for CFHQS z ; 



! 6.4 QUASARS 





CFHQS J0210-0456 


CFHQS J2329-0301 


M BH 


(8.0t^) x 1O 7 M a 


(2.5^1) x 1q8m a 


«MgII 


6.438 ±0.004 a 


6.417 ± 0.002 a 


2[CII] 


6.4323 ±0.0005 




FWHM [cn] 


189 ± lSkms" 1 




I[C1I] 


0.269 ± 0.037 Jykms" 1 


< 0.10 Jykms- 1 b 


L [CII] 


(3.01 ±0.41) x 10 8 L Q 


< 1.1 x 10 s L b 


^00(2-1) 


< 0.014 Jykms" 1 c 




^CO(2-l) 


< 1.6 X 1O 7 L c 




/1.2mm 


120±35MJy 


< 90 /LtJy d 


Lfir 


(2.60 ±0.76) x 10 U L Q 


< 1.9 x 10 11 L d 


SFR 


48 ± 14 M vr" 1 


< 4OM yr- 1 d 


i [CII]/i'FIR 


(1.15 ±0.32) x 10~ 3 




-^CO(1-0)/-^FIR. 


< 8.5 x 10~ 6 





Notes. — 

a Derived from Mgn A2799 observations (Willott et al. 2010b). 
b 3cr upper limit assuming spatially unresolved and line width 
FWHM=300kms- 1 . 

c 3cr upper limit from observations in Wang et al. (2011b) assum- 
ing spatially unresolved and FWZI=300kms~ 1 . 
d 3ct upper limit assuming spatially unresolved. 

data are quoted in Table [T] 

The [CII] line is primarily produced in photo- 
dissociation regions and is strongly dependent on the in- 
terstellar radiation field (Stacey et al. 1991). The ratio 
L[ch]/£fir has an inverse correlation on the radiation 
field strength and has been widely studied at lower red- 
shift. It has been found that the ratio has an inverse cor- 
relation with Lfir (Luhman et al. 2003). Gracia-Carpio 
et al. (2011) found this inverse correlation is even tighter 
if one normalises the far-IR luminosity by the molecular 
gas mass Mh 2 - They attributed this to Lfir/M# 2 be- 
ing more closely related to the physical properties of the 
clouds such as density and temperature. 

Figure m\ shows L^hj/Lfir against Lfir for a com- 
pilation of low (z < 0.4) and high (1 < z < 5) red- 
shift galaxies (Gracia-Carpio et al. 2011 and in prep.). 
The horizontal offset between the low and high redshift 
sources was attributed by Gracia-Carpio et al. (2011) to 
the higher molecular gas mass (for a given Lfir) at high 
redshift (e.g. Tacconi et al. 2010). Also plotted on Fig- 
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Fig. 6. — Ratio of [Cll] to far-infrared luminosity versus far- 
infrared luminosity. Open circles show low-redshift (z < 0.4) and 
filled circles high-redshift (1 < z < 5) galaxies from the collec- 
tion by Garcia-Carpio et al. (2011 and in prep.). Detections of 
quasars at z > 6 are shown individually: J0210-0456 (this paper), 
SDSS J1148+5251 (Maoilino et al. 2005) and ULAS J1120+0641 
(Venemans et al. 2012). Note the outstanding high-redshift Z/fir 
sensitivity of ALMA in this Early Science observation with a mod- 
est integration time. 



ure [6] are data for z > 6 quasars. SDSS J1148+5251 
has a very high Lfir and falls along the sequence of 
2 < z < 5 high-redshift galaxies (Maiolino et al. 
2005). ULAS J 11 20+0641 (Venemans et al. 2012) and 
CFHQS J0210-0456 have more moderate L Fm and fall 
within the region occupied by low-redshift galaxies. In 
the interpretation of the offset being due to higher molec- 
ular gas mass at high-redshift, this would suggest that 
not all high-redshift quasars exist in star-forming hosts 
with higher gas masses than at low redshift. The hor- 
izontal offset that is so striking in Figure [6] is at least 
partially due to a selection effect where previous facili- 
ties did not have the sensitivity to detect more moderate 
Lfir at high redshifts and only ultraluminous contin- 
uum sources were followed up with [C n] observations. It 
is likely there is a large population of hitherto undetected 
high-redshift galaxies with properties like J0210-0456. 

Wang et al. (2011b) reported Very Large Array ob- 
servations aimed at detecting the CO (2—1) emission 
from J0210-0456. The object was not detected and a 
line flux upper limit assuming a full-width zero-intensity 
of 800 kms -1 was reported. We have recalculated the 
line flux limit for the same width as the observed [Cll] 
line (FWZI = 300kms _1 ). The 3cr upper limit on the 
line flux is then < 0.014 Jy km s -1 giving a CO (2 — 1) 
line luminosity limit of £co(2-i) < 1-6 x 1O 7 L0. To 
compare with other works that usually quote the ground- 
state CO transition we assume a luminosity ratio of CO 
(2-1) / CO (1-0) = 7.2 (Stacey et al. 2010; Pa- 
padopoulos et al. 2012). Therefore the CO (1 - 0) limit 
for J0210-0456 is L O O(i-0) < 2 - 2 x 1q6 L and the ratio 
ico(i-o)/-^FiR < 8.5 x 10~ 6 . The limit on this ratio is 
an order of magnitude higher than the values for typ- 
ical ultraluminous high-redshift galaxies and AGN (De 
Breuck et al. 2011) showing that much deeper observa- 
tions are required to detect the molecular gas in galaxies 
such as these, highlighting how [Cn] observations with 
ALMA are the best way to probe the obscured interstel- 
lar medium in typical high-redshift galaxies. 
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Fig. 7. — The background shows the continuum map of J0210- 
0456 and ranges from —3a (purple) to +3<r (yellow). The blue and 
red contours near the central dust continuum source show [C II] 
line emission maps for the blue (—136 < v < — 45kms~ 1 ) and red 
(+45 < v < +136kms- 1 ) wings, respectively. There is an offset 
of 0'.'5 between these peaks, either side of the continuum centroid. 
The inset panel shows a [C II] peak velocity map for those pixels 
with sufficient flux to enable a Gaussian to be fitted. This map 
reveals a smooth velocity gradient across the source. 

3.3. [Cll] dynamics and spatial extent 

Wang et al. (2010) showed that CO line widths of 
z rs 6 SDSS quasars indicate ratios of black hole to galaxy 
mass typically a factor of 10 greater than in the local 
universe. This fits with the results presented in Section 
|3.1| where it was found that z ~ 6 quasar host galaxies 
are growing their black holes at a rate about 10 x faster 
than their stellar mass compared to the local ratio. In 
this work we have measured the [C n] line width for just 
one z = 6.4 quasar, so will only briefly discuss the ratio 
of black hole to dynamical mass at z ~ 6 and defer a 
fuller investigation until more ~ 10 s M Q black hole host 
galaxies have suitable mm interferometry data. 

CFHQS J0210-0456 has a [Cn] line FWHM of 189 ± 
18kms _1 , equivalent to a = 80 ± 8kms _1 for a Gaus- 
sian and ignoring any inclination correction. Based on 
the local relationship (Gultekin et al. 2009) a galaxy 
with a — 80kms _1 would be expected to have Mbh ~ 
3 x 10 6 M Q , a factor of 25 lower than the measured 
-Mbh — 8 x 1O 7 M . This difference is comparable to 
the factor of 10 found by Wang et al. for more luminous 
(Mbh ~ 1O 9 M ) quasars. Even after taking account 
of possible inclination effects (Carilli & Wang 2006) this 
shows the black holes in z 6 quasars to be considerably 
more massive than expected from the local relationship 
between black hole and galaxy mass. 

To fully determine the gas kinematics requires spatially 
resolving the line emission. This should be possible with 
the full ALMA array that will offer spatial resolution 
reaching 20 milliarcseconds (mas) , which is equivalent to 
120 pc. Walter et al. (2009) found that the [C n] emission 
in SDSS J1148+5251 is concentrated within a radius of 
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only 750 pc of the nucleus. This is consistent with the 
small sizes of luminous z s» 6.5 Lyman Break galaxies 
that have effective radii s» 800 pc (Ono et al. 2013). 

The [Cn] line image of J0210-0456 shown in Figure 
[3] shows elongation along an E-W direction, roughly 
aligned with the direction of the beam. The measured 
size of the source is 879 ± 55 mas x 642 ± 75 mas com- 
pared to a beam size of 770 mas x 520 mas. We used 
the CASA IMFIT task to fit a deconvolved model im- 
age to the data. This results in a deconvolved source of 
521 ± 248 mas x 224 ± 297 mas oriented at a position an- 
gle of 128 degrees east of north. The intrinsic source size 
is not well constrained as it is only marginally resolved 
in this image comprised of all fifteen spectral channels of 
the [C11] line. 

A different story emerges when one considers the blue 
and red sides of the [Cll] line separately. Maps were 
made using only the red and blue wings (5 channels each) 
and excluding the centre of the line. Figure [7] shows the 
continuum dust emission as the background image. Su- 
perimposed on this are separate contours for the blue and 
red sides of the [C 11] line. It can be seen that there is a 
spatial offset of 0'.'5 (3 kpc) between these peaks along a 
position angle similar to the 128 deg major axis of the 
marginally resolved full channel image. The inset panel 
of Figure ff\ shows a velocity-centroid map of the [C 11] line 
for pixels with sufficient flux to enable a Gaussian emis- 
sion line to be fitted. There is a clear velocity gradient 
across the source along this same axis with magnitude 
« lOOkrns -1 across a size scale of l'/O (6 kpc). Whether 
this is due to rotation of a galaxy wide disk or has a more 
complex origin in merging multiple components will re- 
quire higher spatial resolution observations. This is simi- 
lar to the [C 11] velocity gradients over this scale observed 
by Wang et al. (2013) in some z ~ 6 SDSS quasars and 
quite different to the compact, intense, central starburst 
observed in SDSS J1148+5251 (Walter et al. 2009). 

4. DISCUSSION 

These observations with ALMA break new ground in 
their sensitivity to moderately star-forming galaxies at 
high- redshift. Even with this sensitivity, only one of 
the two z 6 quasars was detected in line and con- 
tinuum emission. The far-IR luminosity of J0210-0456 
is (2.60 ± 0.76) x 10 n L Q , whereas J2329-0301 remains 
undetected with Lfir < 1.9 x 10 11 Lq, significantly be- 
low the typical far-IR luminosity for a quasar of this 
bolometric luminosity at any redshift. These low far-IR 
luminosities are surprising and place strong constraints 
on the star formation rates in the host galaxies of these 
Eddington-limited quasars. 

In the simplest black hole/galaxy co-evolution sce- 
nario, cosmic stellar mass and black hole mass increase 
in lockstep, ending up at the ratio observed in the local 
universe of M B}i /M stenal = 0.002 (Tundo et al. 2007). 
Detailed simulations show that in individual galaxies the 
phases of star formation and black hole accretion are not 
synchronized (Li et al. 2007) and this accounts for the 
significant scatter of points in Figure[2| It is trivial to cal- 
culate the relationship between star formation rate (lin- 
early related to Lfir) and black hole accretion rate (lin- 
early related to LboI and assuming an accretion efficiency 
of 10%) necessary to achieve Mbh /-^stellar = 0.002. This 
curve is plotted as the dashed line in the upper- left of Fig- 
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ure [2] All the optically-selected quasars observed at mm 
wavelengths are growing their black holes at a relatively 
faster pace than their stellar mass and this is not too 
surprising given that they were selected by their quasar 
emission. Galaxies should lie on both sides of the dashed 
line during their lifetimes in order to reach the local ra- 
tio at z = 0. At z « 2, mm-selected galaxies mostly lie 
to the upper left of the line showing that they are grow- 
ing their stellar mass more rapidly than their black holes 
(Alexander et al. 2005). Lutz et al. (2010) also showed 
that low-luminosity AGN from deep X-ray surveys are 
found on the left side of such a line. 

J2329-0301 is found to be growing its black hole at 
a rate of > 100 x faster than its stellar mass in order 
to reach the local ratio. It has a black hole accretion 
rate of M BH ~ 7M Q yr~ 1 and SFR< 40M Q yr- 1 (3er 
limit, assuming Td < 47 K, no AGN-heated cool dust and 
Salpeter IMF). Khandai et al. (2012) report the results of 
hydrodynamic simulations of z ~ 6 quasar host galaxies. 
These show SFR that range from 100 to 1000 M yr" 1 as 
they evolve during the main black hole accretion phase. 
The simulations are designed to match the properties of 
the most luminous quasars from the SDSS. J2329-0301 
has a black hole accretion rate ~ 3x lower than typi- 
cal SDSS quasars and therefore scaling down the lowest 
simulated SFR by this amount would result in approxi- 
mately the SFR limit observed for J2329-0301. This sug- 
gests that J2329-0301 is observed at a rare phase where 
it has very low SFR compared to its black hole accre- 
tion rate. In the Khandai et al. (2012) simulations, the 
SFR usually drops at the epoch of peak quasar accretion 
due to feedback heating the host galaxy gas. J2329-0301 
appears to have very effectively shut off star formation. 
Although the far-IR data of high luminosity PG quasars 
has many non-detections and the H- ATLAS data are just 
stacked averages (Figure[2]), it would appear that few low 
redshift quasars have such a low ratio of SFR to black 
hole accretion as J2329-0301. This quasar is known to 
be surrounded by a luminous Lya halo at least 15 kpc 
across (Goto et al. 2009; Willott et al. 2011), which sig- 
nifies a huge reservoir of diffuse gas likely photo-ionized 
by the quasar. When this gas cools and accretes on to the 
galaxy, a further bout of star formation is likely. Hayes 
et al. (2013) noted that the Lya emission from galaxies 
with low dust content tends to be more extended than 
that in dusty galaxies, but do not provide a simple ex- 
planation for why this happens. J2329-0301 certainly 
fits this pattern with a low dust content as measured by 
thermal dust emission and a very extended Lya halo. 

The [Cn] line detection in J0210-0456 is narrow 
(FWHM = 189±18 kms" 1 ) and shows only a small veloc- 
ity gradient (« 100 kms -1 ) across a scale of 6 kpc. The 
inclination angle of the [C 11] emission is unconstrained, 
but the narrow line suggests the dynamical mass of the 
system is much lower than would be found in the local 
universe for a galaxy hosting a 10 8 M Q black hole. This 
is in agreement with the results for more massive black 
holes at z w 6 (Wang et al. 2010) and fits with the dis- 
cussion above that shows optically-selected z rs 6 quasars 
have experienced enhanced black hole accretion relative 
to their stellar mass accumulation. 

These early observations with ALMA are a prelude to 
increased resolution and sensitivity observations to come. 
A larger sample of z sa 6 quasars with a wide range of 
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black hole masses is necessary to get a more complete 
picture of the relationship between black hole growth and 
star formation. Spatially resolved observations on scales 
of ~ 100 pc will reveal the dynamical state of the star 
forming gas and enable more accurate determination of 
dynamical masses. 

Thanks to Javier Gracia-Carpio for providing unpub- 
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ies. Thanks to staff at the North America ALMA Re- 
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